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We consider a many-body generalization of the Kapitza pendulum: the periodically-driven sine-Gordon 
model. We show that this interacting system is dynamically stable to periodic drives with finite frequency and 
amplitude. This finding is in contrast to the common belief that periodically-driven unbounded interacting sys¬ 
tems should always tend to an absorbing infinite-temperature state. The transition to an unstable absorbing state 
is described by a change in the sign of the kinetic term in the effective Floquet Hamiltonian and controlled by 
the short-wavelength degrees of freedom. We investigate the stability phase diagram through an analytic high- 
frequency expansion, a self-consistent variational approach, and a numeric semiclassical calculations. Classical 
and quantum experiments are proposed to verify the validity of our results. 

PACS numbers: 67.85.-d, 03.75.Kk, 03.75.Lm 


I. INTRODUCTION 

Motivated by advances in ultra-cold atom^HIl, the stabil¬ 
ity of periodically-driven many-body systems is the subject 
of several recent studie^^^^^. According to the second law of 
thermodynamics, isolated equilibrium systems can only in¬ 
crease their energy when undergoing a cyclic process. Eor 
many-body interacting ergodic systems, it is often assumed 
that they will heat monotono usly, a symptotically approaching 
an infinite-temperature statJ^KSEH jn contrast, for small sys¬ 
tems such as a single two-level system (spin), thermalization 
is not expected to occur and periodic alternations of heating 
and cooling (Rabi oscillations) are predicted. A harmonic os¬ 
cillator can display a transition between these two behaviors, 
known as “parametric resonance’® depending on the ampli¬ 
tude and frequency of the periodic drive, the oscillation am¬ 
plitude either increases indefinitely, or displays periodic oscil¬ 
lations. An interesting question regards how much of this rich 
dynamics remains when many-degrees of freedom are consid¬ 
ered. 

This question was addressed for example by Russomanno 
et alW, who studied the time evolution of the transverse-field 
Ising (Tl) model. This model is integrable and never flows 
to an infinite-temperature state. In the translational-invariant 
case, this result can be rationalized by noting that the Tl model 
is integrable and can be mapped to an ensemble of decoupled 
two-level systems (spin-waves with a well defined wavevec- 
tor), each of whom periodically oscillates in time and never 
equilibrates. In this sense, the findings of Ref.s 09I11I12I14I 
on periodically-driven disordered systems subject to a local 
driving falls into the same category: many-body localized 
(MBL) systems are effectively integrable because they c an be 
described as decoupled local degrees of freedom as welP®. 
When the driving frequency is higher than a given threshold. 


MBL systems remain localized and no not thermalize. This 
is in contrast to ergodic systems, which are expected tother- 
malize to an infinite temperature for any periodic driv^^®3111 
These find ings are in apparent contradiction to earlier numer¬ 
ical studiePiSlEI] found indications of a finite stability 
threshold in non-integrable systems as well. 

To investigate this problem in a systematic way, we con¬ 
sider here a many-body analog of the Kapitza pendulum: 
the periodically-driven sine-Gordon model. This model 
is well suited for analytical treatments, including a high- 
frequency expansion, quadratic variational approaches, and 
renormalization-group methods. Unlike previously-studied 
spin systems, the present model has an unbounded single¬ 
particle energy spectrurrP^. Thanks to this property, infinite- 
temperature ensembles are characterized by an infinite energy 
density and are easily identified. We show the emergence of a 
sharp “parametric resonance”, separating the absorbing (infi¬ 
nite temperature) from the non-absorbing (periodic) regimes. 
This transition survives in the thermodynamic limit and leads 
to a non-analytic behavior of the physical observables in the 
long time limit, as a function of the driving strength and/or 
frequency. We conjecture that this transition corresponds to a 
mean-field critical point of the many-body Eloquet Hamilto¬ 
nian. Our finding enriches the understanding of the coherent 
dynamics of parametrically forced system and paves the road 
toward the search of unconventional dynamical behavior of 
closed many-body systems. 


II. REVIEW OE A SINGLE KAPITZA PENDULUM 

Before entering the domain of many-body physics, we 
briefly review the (well understood) case of a classical sin¬ 
gle degree of freedom. We consider the Hamiltonian of a 
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FIG. 1: Stability diagram of the classical Kapitza pendulum, adapted 
from Ref. 1281 . In the colored areas at least one fixed point is stable, 
while in the white areas both minima are unstable and the system 
is fvariational methodsully ergodic. The parameters go, g\ and 7 
are defined in Eq. l[^ as g{t) = 50 + <?i cos( 7 t). The green lines 
correspond to the stability threshold of the first parametric resonance, 
Eq. 0 



periodically-driven simple (non-linear) pendulum, also known 
as Kapitza pendulurrP^l and described by the Hamiltonian 

Hit) = - g{t) cos(^), with g{t) =50+51 cos{-ft) . 

( 1 ) 

Here p and (j) are canonically-conjugated coordinates satisfy¬ 
ing {p, 0 } = —i, where {•,•} are Poisson brackets. For 51 = 0 
the system displays two classical fixed points: a stable one at 
^ = 0 and an unstable one at (p = tt. (Throughout this paper 
we assume without loss of generality that po > 0.). 

In the presence of a periodic drive (51 7^ 0 ), the unsta¬ 
ble fixed point can become dynamically stable. This coun¬ 
terintuitive result was first obtained by Kapitzd^ in the high- 
frequency limit 7^ ^ 50) 5i- By averaging the classical equa¬ 
tions of motion over the fast oscillations of the drive, Kapitza 
found that the “upper” extremum (j> = tt becomes stable for 
large enough driving amplitudes g^ > go^"^ 12 . This pioneer¬ 
ing work initiated the field of vibrational mechanicJ^, and 
the Kapitza’s method is used for description of periodic pro¬ 
cesses in atomic physicJ^, plasma physicJ^, and cybernetical 
physicJ^. 

For finite driving frequencies 7^ ~ po: 5 i the lower fixed 
point {<p = 0 ) can become dynamically unstable as well. This 
transition can be analytically studied for example by applying 
the quadratic approximation cos{cj)) — 1 — 0^/2 to Eq. (IH 
(valid for small pi). The resulting Hamiltonian corresponds 
to a periodically-driven harmonic oscillator, with frequency 
Wo = driving frequency 7, and driving amplitude pi. 

For infinitesimal driving frequencies (pi —>^ 0 ), this system 
displays parametric resonances at 7 = 2 uJo/n = 2 y/go/n, 
where n is an integei^^. For finite driving amplitudes each 
resonance extends to a finite region of driving frequencies® 
in particular the first parametric resonance (n=l) extends to 

251 < - 4 po| • ( 2 ) 

Because the subsequent parametric resonances occur at lower 
frequencies, one finds that the point ^ = 0 is always dynam¬ 
ically stable for driving frequencies that are larger than the 
critical value 7c = V^ 4 po + 2pi, while regions of stability and 


instability alternate for 7 < 7c. In particular, the system is sta¬ 
ble for po/7^ > 0.25 and small driving amplitudes. In what 
follows we will refer to these two stability regions as “large 
driving frequency” and ’Targe po” respectively. As a side re¬ 
mark, we note that this quadratic approximation is valid for 
the quantum case as well, provided that the initial state is close 
to the \(j> = 0 { state. In this case, the resulting stability phase 
diagram is expected to be the same. 

Subsequent numerical studies of the classical Hami ltonian 
([^ lead to the stability diagram reproduced in Fig. [ijpl In this 
plot, the white areas represent regions in the parameter space 
in which both extrema are unstable and the system flows to¬ 
wards an infinite-temperature state, independently on the ini¬ 
tial conditions. In contrast, in the colored regions at least 
one of the two extrema is stable, and the system is not er¬ 
godic. In this case the pendulum can be confined to move 
close to one of the extrema and will not in general reach 
a steady state described by an effective infinite temperature. 
The stability of the lower fixed point (green line in Fig. [T]) is 
well approximated by the boundaries of Eq. In particular, 
for go = 0 (where the upper and lower extrema are math¬ 
ematically equivalent), the system is dynamically stable for 
5 i < 5 c ~ 0.457^, and dynamically unstable for pi > pc. 
As a main result of this paper we will show that this transition 
remains sharply defined even in the many-body case. 


III. MANY-BODY KAPITZA PENDULUM 


To explore the fate of the dynamical instability in a many- 
body condition we consider an infinite number of coupled 
identical Kapitza pendula, depicted in Eig.j^a). This system 
is described by the periodically-driven Erenkel-Kontorovd® 
model 


H n 1 , , , , q(t) , , , 

- Y ^ cosip,) , 

( 3 ) 

where Pi, pi are unitless variables satisfying {Pj, pk] = 
—i( 5 j fc, and p(f) is defined in Eq. Q. The energy scale A de¬ 
termines the relative importance of the coupling between the 
pendula and the forces acting on each individual pendulum: 
in the limit A —^ 0 we expect to recover the case of an iso¬ 
lated periodically-driven pendulum. In the continuum limit, 
the model 0 can be mapped to the periodically-driven sine- 
Gordon model 


iT = AV 


H = dx 


“ 3{t) cos,[p) 


( 4 ) 


P{x) and p{x) are canonically-conjugate fields, 
{P{x), p{x')} = —i 6 {x — x') and K is the Luttinger param¬ 
eter (we work in units for which the sound velocity is m = 1 ). 
The parameter A enters as an ultraviolet cutoff, setting the 
maximal allowed momentum: p{x) — j^f^dq/ipiir) e'^'^^pq. 
The model 0 can also be realized using ultracold atoms 
constrained to cigar-shaped traps. In this case, the field 
p = Pi — p2 represent the phase difference between the two 
condensates, and the time-dependent drive can be introduced 
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FIG. 2: Proposed physical realizations, (a) A one-dimensional ar¬ 
ray of Kapitza pendula, which are coupled and synchronously-driven 
realizes the time-dependent Frenkel-Kontorova model 0- (b) Two 
coupled one-dimensional condensates, whose tunneling amplitude is 
periodically driven in time, realize the time-dependent sine-Gordon 
model 0- More details about the experimental realizations are pre¬ 
sented in Sec. m 


by periodically modulating the transversal confining poten¬ 
tial, as shown in Fig. |^b). Because the helds P and (p are 
continuous variable, the energy densities of the Hamiltonians 
0 and 0 are unbounded from above and allows to easily 
distinguish an absorbing behavior (in which the energy 
density grows indehnitely in time) from a periodic one. This 
situation differs from previously-considered spin models, 
whose energy density is generically bounded from above. 


IV. INFINITE FREQUENCY EXPANSION 

To understand the effect of periodic drives on dynamical in¬ 
stabilities and localization, it is convenient to dehne the effec¬ 
tive Hamiltonian (often termed “Floquet Hamiltonian”) iJeff 
as 

U{T) = , (5) 

where U (T) is the evolution operator over one period of time 
T = 2^!^. In the case of the harmonic oscillator, it has been 
shown that the parametric resonance can be easily understood 
in terms of the effective HamiltoniarP^ in the stable regime 
the eigenmodes of H^s are normalizable and any initial state 
can be expanded in this basis, leading to a periodic dynam¬ 
ics. In contrast, in the unstable regime the eigenmodes of iFeff 
become not normalizable, in analogy to equilibrium Hamilto¬ 
nians that are not bounded from below (such as iF = —p^), 

and the dynamics become s abso rbing. 

The Magnus expansiorP^^ is an analytical tool to de¬ 
rive the effective Hamiltonian in the limit of large drive fre¬ 
quency. The hrst-order term of the Magnus expansion, iFeff = 


(1/r) fj' dt H{t), has a simple physical interpretation: when 
the driving frequency is inhnite, the system perceives only the 
time-averaged value of H{t). This observation has been suc¬ 
cessfully employed, for example, to engineer optical lattices 
with negative tunneling amplitudd^S. In our case, the aver¬ 
age Hamiltonian is simply described by the time-independent 
sine-Gordon model (Eq. 0 with gi = 0 ) and the drive has 
no effect on the system. The third-order Magnus expansion 
delivers: 

f dx 

Heff = Hll+H'+H''- I ^ [poCOs((/))+5COs(2(;ii)] , (6) 

where H is the Luttinger liquid Hamiltonian (Eq. Q with 

9i= 90 = 0 ), 



H" = -g" J cos{(j)); ( 8 ) 


9' = 5 " = ^ and 5 = (|pi - go). See Ap- 

pendixj^for the details of this derivation. Eq. (j^ is analogous 
to the effective Hamiltonian of the single Kapitza pendulurrPl 
The stability of the “upper” extremum p = tt is captured 
by the interplay between go cos{ 4 >) and gcos{ 2 (j)): this point 

becomes dynamically stable when go < Ag, or equivalently 

2 

when go < taking into account that the Magnus expan¬ 
sion is valid in the limiP^of 50 / 7 ^ ^ 1- 

The term H' leads to the dynamical instability of the sys¬ 
tem: having a negative sign, it suppresses the kinetic energy 
~ P^, eventually leading to an inversion of its sign. Using 
a quadratic variational approach, we can approximate H' « 
—4g'f dx{cos{(l>))P^. The stability of the cp = 0 minimum 
can then be studied through {cos{p)) « 1 , leading to a renor- 
malizedkinetic energy iTP^/ 2 —Ar^(pi/ 7 ^)P^ = /% 

with 

iTeff = k{i- (9) 


This expression indicates that the system is dynamically stable 
for AT( 71 / 7 ^ < 0.5. Eor larger driving amplitudes (or smaller 
driving frequencies) the kinetic energy becomes negative and 
the system becomes unstable. 

Higher-order terms in the Magnus expansion can 
be used to determine the qualitative dependence of 
the critical driving amplitude as a function of A /7 
and Kgi/^'^. Eor example, the hfth-order term of 
the Magnus expansion contains terms proportional 
to 51 (f)/ 74 {(a./>) 2 , {(9.0)2, {P^{P^cos</}}}} ~ 


P^d^ cos( 0 ), which renormalizes g' by a factor of 8^ j . 
This positive contribution leads to a decrease of the critical 
driving amplitude as a function of A/ 7 . As mentioned above 
(see Eq.0), A sets the coupling between the Kapitza pendula 
and is indeed expected to shrink the stability region of the 
system. 

Since the Magnus expansion generates an inhnite number 
of terms, one may suspect that the full series could renormal¬ 
ize the critical amplitude to zero (making the system always 
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dynamically unstable). To address this point, we now resort 
to the quantum version of the problem, where powerful ana¬ 
lytical techniques are available. 


V. IfENOIfMALIZATION GROUP ARGUMENTS 

Let us now consider the quantum version of the Hamilto¬ 
nian Q. In the absence of a drive {gi = 0), this model cor¬ 
responds to the celebrated sine-Gordon modePSEl. Its ground 
state displays a quantum phase transition of the Kosterlitz- 
Thouless type: for K > Kc = Sir — o{g/A) the cosine term 
is irrelevant and the model supports gapless excitations, while 
for K < Kc the system has a finite excitation gap A. The 
presence of a gap is known to have significant effects on the 
response of the system to low-frequency modes, by exponen¬ 
tially suppressing the energy absorption. In contrast, at large 
driving frequencies the excitation gap is expected to have little 
effect. 

To analyze the limit of large driving frequencies, we pro¬ 
pose to consider the Floquet Hamiltonian, as defined by the 
Magnus expansioJ^. Specifically, we consider the ground- 
state properties of this Hamiltonian, which can be conve¬ 
niently studied through the renormalization group (RG) ap¬ 
proach. The existence of a well-defined ground state for the 
Floquet Hamiltonian implicitly demonstrates the ergodicity of 
the system: unstable systems are genetically expected to have 
non-normalizable eigenstateJSl 

In our case, the first-order term of the Magnus expansion 
corresponds to the well-known sine-Gordon model. When 
Keff > Stt this Hamiltonian flows under RG towards the Lut- 
tinger liquid theory (Eq.Q with g{t) = 0). Higher order 
terms are given by commutators of the Hamiltonian at differ¬ 
ent times (see Appendix [A|). Because the time-dependent part 
of our Hamiltonian is proportional to cos{(j)), each commuta¬ 
tor necessarily includes cos{(j)), or its derivatives. With respect 
to the Luttinger-liquid fixed point, these terms are irrelevant in 
an RG sense, and are not expected to affect the ground-state- 
properties of the Floquet Hamiltonian. If this is indeed the 
case, the (asymptotic) expectation values of physical operators 
such as 1P and |Pg^ are finite (and proportional to Kgg/\q\ 
and \q\/KcS respectively), indicating that the system does not 
always flow to an infinite-temperature ensemble. When the 
frequency is reduced, the amplitude of higher-order terms of 
the Magnus expansion increases: although irrelevant in an RG 
sense, if sufficiently large, these terms can lead to a transition 
towards an unstable regime. 

Many-body quantum fluctuations have an important effect 
on the stability of the inverted pendulum (/) = tt as well. As 
mentioned above, this effect is related to the interplay be¬ 
tween cos{(j)) and cos(20) in Eq. In the ground-state of 
this Hamiltonian, quantum fluctuations change the scaling di¬ 
mension of an operator cos{acj)) to 2 — a^K/A. Thus for fi¬ 
nite Kes, the term with g is less relevant than the term with 
go from a renormalization group (RG) point of view, making 
the upper extremum less stable than in the case of a simple 
pendulum. A simple scaling analysis reveals that the stabil¬ 
ity boundary between the two extrema at 0 and tt is renor¬ 


malized by approximately (g/A)^, where A is the theory cut¬ 
off. This is a strong indication that, for any finite A, the ex¬ 
tremum at TT can still be made stable with large enough driving 

amplitudeJSl. 

VI. QUADRATIC EXPANSION 

In analogy to the single Kapitza pendulum, the simplest 
way to tackle the many-body Hamiltonian Q is to expand 
the cosine term to quadratic order. In particular, to analyze 
the stability of the ^ = 0 minimum we can use cos{(j)) —>■ 
1 — (j)'^/2. In this approximation the system becomes equiva¬ 
lent to a set of decoupled harmonic oscillators with Hamilto¬ 
nian H = Hq, where 

Hq = ^Pq + ^ ( 9 ^ + -^50 + Kgi cos(7f)) . (10) 

Eq. ( [T0) l corresponds to the Hamiltonian of a periodically 
driven harmonic oscillator, whose stability diagram is well 
known. In the limit of gi 0, the system is dynamically 
stable as long as 7 > 2max[a;q] = 2 max[A/g2TiI~/^^] = 
2-\/A2 -I- Kgo, or 7 < 2min[wg] = 2min[i/q2 _|_ Kgo] = 
2y/KgQ. In analogy to Eq.(j^, for finite gi the stability condi¬ 
tion is modified to 

2Kgi < max [ 7 ^ - 4 (Kgo -f A^) , 4 (Kgo) - 7 ^] (I I) 

The resulting dynamical phase diagram is plotted in Eig. 

In the limit of A 0 we recover the stability diagram of 
the lower minimum of an isolated Kapitza pendulum: this is 
demonstrated by the quantitative agreement between the green 
curves of Eigj^and Eig. 


00 


FIG. 3: Stability diagram of the periodically-driven sine-Gordon 
model, as obtained from the quadratic expansion of the lower min¬ 
imum. The system displays two distinct stability regions at large 7 
and large go, respectively (see text). The green line refers to the limit 
A/7 —>■ 0, where the problem maps to the stability of the lower min¬ 
imum of the Kapitza pendulum (green line of Fig.[TJ. 





















5 


For finite A we observe two distinct stability regions, char¬ 
acterized by a different dependence on the ultraviolet cut¬ 
off A. The first stability region is adiabatically connected 
to the limit of infinite driving-frequency, 7 —oo (the ori¬ 
gin in Fig. [^). Its boundaries are described by -b 

^ 50 / 7 ^ + ^ = 0.25. For po —^ 0 and A —>• 0 this ex¬ 
pression precisely coincides with the result obtained from the 
large-frequency Magnus expansion, Kgij^'^ = 0.5. The sta¬ 
bility region is strongly suppressed by A and eventually dis¬ 
appears at A /7 « 0.5. This indicates that it is related to the 
stability of the degrees of freedom at the shortest length scales 
A“^: its character is therefore predicted to be analogous to the 
stability of a single Kapitza pendulum. In contrast, the second 
stability region, at large g^ is roughly independent on A/ 7 . 
This stability region is determined by the resonant excitation 
of the lowest frequency collective excitation of the system in 
the static (7 —> 0 ) regime, where the system is dynamically 
stable due to the presence of a gap, A. In the quadratic ap¬ 
proximation, the gap is given by A « y/Kg^, and the stability 
condition 7 < 2A reads Kg^/^"^ < 0.25. 

The present quadratic expansion bares a close resemblance 
to the analysis of Pielaw^, who considered periodic modula¬ 
tions of the sound velocity (but keeping po = 0). The validity 
of these quadratic approximations is however undermined by 
the role of non-linear terms in non-equilibrium situations (and 
specifically noisy environment^ and quantum quencheJSl). 
Even if irrelevant at equilibrium, non-linear terms enable 
the transfer of energy between modes with different momen¬ 
tum and are therefore necessary to describe thermalization. 
The question which we will address here is whether mode¬ 
coupling effects are sufficient to destroy the dynamical insta¬ 
bility described above. 


VII. SELF-CONSISTENT VARIATIONAL APPROACH 


The above-mentioned quadratic approximation can be im¬ 
proved by considering a generic time-dependent Gaussian 
wavefunction, along the lines of Jackiw and KermarP5 


'I'„[(/(x)] = Tlexp 






( 12 ) 


Eq. ( |T 2 ] i is a particular case of the generic wavefunction pro¬ 
posed by Cooper et a/.EH, valid when the expectation value 
of the field and its conjugate momenta are zero: {4>{x)) = 
{P{x)) = 0. The operator A ~ (detG)^/^ is the normal¬ 
ization constant to ensure the unitarity of the evolution at all 
times: 


= J 'i'vm = 1 . ( 13 ) 


The functions Gx,y and ^ are variational parameter to 
be determined self-consistently. To this end, we invoke the 
Dirac-Erenkel variational principle and define an effective 
classical Lagrangian density as follows: 


Eeff = J !?[</] [idt - H[<P,d/d4>]) (14) 


Eor a translation-invariant system, we find Gx,y = Gx-y = 
l/(27r) dk Gk where A is the UV cutoff, and the 

effective action SeS = f f dx £eS is given by 


Sci = Z(t) g(t) + 



where 


-^KG^^- 2 KSkGkSk - ^ Gfe), 


(15) 


The equations of motion are given by the saddle point of 
the effective action. By requiring SSen/SG = 5Ses/5Yj = 0 
we obtain 


Gfc=4iTGfcSfe, (16a) 

tk = \KGl^-2KYl-^-\z{t)g{t) (16b) 


In the following calculations we assume the system to be 
initially found in a stationary state satisfying 


Gk = ^ 


1 


2 

where Ag is self-consistently given as: 


Ag = goK exp - 


1 


dkK 


2 7 _a Stt 2 _|_ ^2 j 


(17) 


(18) 


Assuming Ag A, the above equation gives Ag « 
(( 7 oA'/ 2 )[Ag/( 2 A)]^/^'^, which implies a critical point at 
Kc = Stt (Kosterlitz-Thouless transition). The cosine is rele¬ 
vant (irrelevant) for K < Kc {K > Kc). 

In this initial gapped phase, the classical oscillation fre¬ 
quency, y/Kgo is renormalized by the factor Z due to quan¬ 
tum fluctuations (see Eqs. and ([T 6 |)). The introduction 
of a modulation to the amplitude of the bare cosine poten¬ 
tial leads to one of the two following scenarios (i) Unstable 
(ergodic) regime: the driving field amplifies quantum fluctua¬ 
tions (i.e. leads to “particle generation” via parametric reso¬ 
nance) and closes the gap, i.e. Z{t) 0. Once the gap closes, 
it remains closed; we take this as an indication of the run¬ 
away to the infinite-temperature limit (we can also study the 
absorbed kinetic energy in this formalism as well); (ii) Stable 
(non-ergodic) regime: quantum fluctuations remain bounded, 
Z{t) stays finite at all times, and (p remains localized. These 
two regimes are indicated in Eig. as white and steel-blue 
regions, and are in quantitative agreement with the results of 
the quadratic approximation, Eig. The apparent inconsis¬ 
tency for small A /7 <C 1 and Kgijj^ <C 1 is simply due to 
finite-time effects (see also Appendix [B]l 

The effects of quantum fluctuations is analyzed in Eig.Qb). 
This plot displays the stability diagram for a fixed and small pg 
{goKj^f^ = 10“^) as a function of pi, AT and A. Eor small K 
we reproduce the large-frequency stability lobe of Eig. |^a). 
Indeed, the quadratic approximation is expected to become 










6 


exact in the limit of X —0. Finite values of the Luttinger pa¬ 
rameter K significantly shrink the volume of the non-ergodic 
(stable) regime. This result indicates that many-body quan¬ 
tum fluctuations promote ergodicity. At the same time, our 
analysis indicates that a finite region of stability can survive 
in the thermodynamic limit, even in the presence of quantum 
fluctuations. 


(a) K = O.Itt 



VIII. SEMICLASSICAL DYNAMICS 


To further demonstrate the existence of an instability tran¬ 
sition at a finite driving frequency, we now numerically solve 
the classical equations of motion associated with the Hamil¬ 
tonian 0 - Specifically, we focus here on the stability re¬ 
gion at large driving frequencies and along po = 0 (magenta 
curve of Fi g.s and Qa)). Following the truncated-Wigner 
prescriptiorPSEsT we randomly select the initial conditions 
from a Gaussian ensemble corresponding to the ground state 
of 0 with g{t) = 0 and solve the classical equations of mo¬ 
tion associated with the Hamiltonian Q. Although not shown 
here, we checked that other choices of initial conditions lead 
to qualitatively similar results. Fig.|^shows (a) the time evo¬ 
lution of the average kinetic energy and (b) of the expectation 
value {cos{(j>)) for different driving amplitudes. The former 
displays a sharp increase in correspondence of the expected 
dynamical transition. Fig. |^a) shows the average kinetic en¬ 
ergy Eao, and the oscillation amplitude 6 cos{(j)) as a function 
of driving frequency: Although is smooth as a function 
of the driving amplitude, its first derivative (Fig.j^b)) shows a 
sharp kink at a critical value of iTpi/ 7 ^. Non-discontinuities 
in the second derivative of the energy are clear evidence of 
continuous transitions. In contrast, Scos{(j>) presents a kink 
itself at the critical value of Kgi/j'^. 

From the position of the kink in either Fig. |^b) or (c) we 
compute the dynamical phase diagram shown in Fig. We 
And that larger K lead to a reduction of the stability region. 
A similar result was obtained using the self-consistent varia- 


(b) A50/7" = 10"'^ 



2 ' 0.6 


(a) 



0 10 20 30 

t {y/2K} 

(b) 



0 10 20 30 

t (y/2n) 


-Kg//=0.62 

- Kg,//=0.47 
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- Kgy/=0.11 

- Kg^//=0.046 
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FIG. 4: The dynamical regimes (a) as a function of 

{goKly'^, giKjy'^, Af^) for fixed K = 0.1, (b) as a function 
of {K, giK/y^, A/y) for fixed goKj^'^ = 10“"*. Blue (white) 
regions correspond to points in the parameter space where the sys¬ 
tem is stable (unstable). Solid lines in (a) correspond to the results 
of the quadratic approximation (Fig. and identify two distinct 
stability regions, respectively at large 7 and large go. In this plot 
the stability criterion is arbitrarily set to Z{Tf)jZ{0) > 0.95 with 
Tf = 100(27r/7) (See also Fig. |^for details about the finite-time 
scaling.) 


FIG. 5: (a) Time evolution of the normalized kinetic energy 

O'kin(f) = Ekin(f)/i5kin(f = 0 ) - 1 , with i?km(t) = 

for K = O.dvr, go = 0, A /7 = 0.04, L = 200, 
N — 400. For small drives and large frequencies (lower curves) 
the system is stable and periodically oscillates with a period tt/A = 
12.5(27r/7). Upon reaching a critical value of A'gri/ 7 ^ the oscil¬ 
lations are substituted by an exponential increase of the energy. The 
dashed line is a guide for the eye. (b) Time evolution of (cos(2(/)) for 
the same parameters as before. The oscillations become very large 
around a critical value of Kgxj'y^. 
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(c) 

0.151-^^- 




Kg/Y^ 



- A/-lt0.032 

- A/-lt0.052 

- A/-lt0.085 

A/-lt0.14 

- A/-it0.23 

A/y=037 


FIG. 6: (a) Normalized kinetic energy (Jkin at long times (t = 
100 X 2-kI'^') as a function of the normalized driving amplitude 
for K = O.drr, go = 0, 1/ = 200. This quantity displays 
a sharp kink around a critical value indicated hy the dashed line, (b) 
First derivative of the kinetic energy, showing a sharp discontinu¬ 
ity in its derivative, confirming the hypothesis of a second-order-like 
phase transition, (b) Oscillation amplitude of {cos(2<^)) at long times 
(T = 30 X 27r/7) as a function of the normalized driving amplitude 
Jfgi / 7 ^. This quantity displays a sharp peak around a critical value, 
corresponding to the sharp increase in the kinetic energy identified in 
(a). With increasing A /7 the peak moves to lower values of the drive 
amplitude and becomes less pronounced. 



FIG. 7: Stability diagram; critical drive as a function of the UV cut¬ 
off A /7 for go = 0, as obtained through the present semiclassical 
approximation (-I-). The dashed line is a guide for the eyes. In the 
limit A /7 —>■ 0 all the curves tend towards the critical value of a 
single Kapitza pendulum, A'gi/ 7 ^ ~ 0.45. The error bars refer to 
the numerical uncertainty in the position of the peak and demonstrate 
that the peaks do not broaden and remain well defined for finite A/ 7 . 
The magenta curve corresponds to the result of the quadratic expan¬ 
sion (Eq. and magenta curve in Fig.[^. 


tional approach (Fig. |^. In particular, the stability diagram 
is mainly insensitive to the initial conditions. Specifically, we 
repeated the semiclassical dynamics for a different set of ini¬ 
tial conditions (corresponding to an initially gapped state with 
local correlations only) and observed a qualitatively similar 
dynamical phase diagram. 

The above calculations refer to physical observables mea¬ 
sured after a finite number of driving periods. One important 
question regards the evolution of the phase diagram of Fig. 
as a function of time. In particular, one may wonder whether 
the stability region gradually shrinks and disappears in the 
infinite-time limit. In other words, can the dynamical tran¬ 
sition be induced by applying an infinitesimal drive for very 
long times? To address this question, in Fig. we plot the 
critical driving amplitude as a function of the number of driv¬ 
ing periods. Although the critical driving amplitude initially 
decreases as a function of time, we observe that at long times 
it tends to a finite asymptotic value; the stable regime occu¬ 
pies a finite region in the parameter space even in the asymp¬ 
totic long-time limit. The number of driving periods plays an 
analogous role to the size of the system in equilibrium phase 
transitions: when appropriately rescaled, any physical quan¬ 
tity is associated with a well-defined asymptotic scaling limit 
(see Appendix 1^1. 

The present calculations demonstrate the existence of a lo¬ 
calized (non-ergodic) phase in the thermodynamic limit. Fol¬ 
lowing D’Alessio et alV\ this result can be interpreted as a 
(many-body) energy localization transition. Remarkably, the 
present semiclassical approach involves the solution of clas¬ 
sical equations of motion (subject to quantum initial condi- 



FIG. 8 : Critical value of the driving field as a function of the waiting 
time Tfln. The critical field is determined by two independent meth¬ 
ods, namely by observing the peaks of (jk as in Fig. |^b), and of 
5 cos, as in Fig. |^c). Both methods show an initial decrease of the 
critical drive, followed by a saturation at a finite asymptotic value. 
Numerical values: go = 0, A = O.Ttt, A /7 = 0.1 


tions). Our findings are in contrast to an earlier conjecture 
formulated in the context of MBL systems with randomness 
by Oganesyan et aZ.® and suggesting that such a transition has 
a pure quantum origin and does not occur in classical systems. 


IX. EXPERIMENTAL REALIZATION AND OUTLOOK 

We now offer more details about the two physical realiza¬ 
tions that were anticipated in Fig.|^ The two proposed expert- 
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ments belong respectively to the classical and quantum realm. 
The present calculations suggest that the two model should 
have a qualitatively similar stability diagram, but an experi¬ 
mental verihcation is required. The realization of the Hamilto¬ 
nian 0 using classical elements seems particularly appealing 
and relatively easy: this model describes an array of pendula 
attached to a common periodically oscillating support, and 
coupled through nearest neighbor couplings (see Fig|^a)). 

A quantum version of the periodically-driven sine-Gordon 
model (|^ can be realized using ultracold atom^^oHSl gy 
trapping the atoms in cigar-shaped potentials it is possible 
to obtain systems in which the dynamics is effectively one 
dimensional. Specifically, the transversal conhnement can 
be generated through laser standing waves (see for example 
Ref 154)), or through magnetic helds induced by currents 
running on a nearby chip (see for example Ref. ED). In 
both cases, the amplitude of the transverse confinement can 
be easily modulated over time, allowing to realize the setup 
shown in Fig.j^b). Here a time-dependent transverse confine¬ 
ment induces a time-dependent tunneling coupling between 
two parallel tubes. In model Q this coupling is modeled 
as f dxcos{cj)), where (j){x) = (f>i{x) — (p^j x) is the local 
phase difference between the quasi-condensatesSEHESl Deal¬ 
ing with a system of two coupled tubes, the Luttinger param¬ 
eter is doubled with respect to a single tube and equals to 
K = 2>/r, where F = mA/fi^pQ, A = p. is the chemical 
potential, m is the mass of the atoms, and po their average 
one-dimensional densit}^^. Atoms on a chip are characterized 
by relatively small interaction energies, constraining the max¬ 
imal value of the achievable K. Typical experimental values 
of 7 are in the order of 7 < 10“^, or AT < 0.2. The parame¬ 
ter g{t) is set by the instantaneous (single-particle) tunneling 
rate through Jj^ [t) = Kg{t) /p. These experiments are there¬ 
fore constrained to g(t) > 0, or gi < g^. Fig.s|^ andshow 
that the a wide region of stability is expected in the physically 
relevant regime (small K and gi < go), demonstrating the 
feasibility of the proposed experiment. 

An alternative procedure to periodically drive the two cou¬ 
pled tubes involves the modulation of the chemical poten¬ 
tial difference between them 6p{t) = pi(f) — P 2 (f) (while 
keeping approximately fixed the tunneling element). In the 
bosonization language of Eq.Q, this corresponds to a time- 
dependent field that couples to the atomic density. An appro¬ 
priate gauge transformation allows one to map this problem 
into a phase-modulated sine-Gordon model in which the drive 
enters through the phase drive as gcos{(p — 4>o{t)). We 
postpone the detailed analysis of this model to a future publi¬ 
cation: preliminary calculations indicate that the stability di¬ 
agram is analogous to Fig. However, we do not know yet 
whether the two models belong to the same universality class. 

The sine-Gordon model is often discussed in the context 
of isolated tubes under the effect of a longitudinal standing 
waves, or optical lattice, as well. At equilibrium this model 
describes the Luttinger liquid to Mott insulator quantum phase 
transition of one dimensional degenerate gases (see Ref.for 
a review). Interestingly, the inverted pendulum ^ = tt cor¬ 
responds to a distinct topological phase, the Haldane insula- 
toiE 2 E 2 ] At equilibrium non-local interactions are needed in 


order to stabilize this phase. As explained above, in the pres¬ 
ence of a periodic drive it may be possible to stabilize it by 
tuning the system in the region where only the “upper” ex¬ 
tremum is stable. This situation is analogous to the recent 
proposal of Greschner et al^. In relation to the dynamical 
stability of this phase, one however needs to notice that the 
optical lattice affects the Luttinger parameter as well. A cor¬ 
rect description of these experiments therefore involves the 
additional transformation K —>■ K{t) = Kq + Ki cos{'jt), 
with possible signihcant consequences for the resulting sta¬ 
bility diagram. Finally, it is worth mentioning other possible 
realizations of the model Q, including for example RF cou¬ 
pled spinor condensateand arrays of Josephson junctionJ^ 
in the presence of time-dependent magnetic fields . 

To summarize, in this paper we combined different analyt¬ 
ical and numerical tools to study the periodically-driven sine- 
Gordon model. We applied a controlled high-frequency ex¬ 
pansion, the Magnus expansion, to derive an effective (Flo- 
quet) Hamiltonian. Employing ideas from the renormaliza¬ 
tion group (RG) method, we propose the existence of a non¬ 
absorbing (non-ergodic) fixed point, in which the system is 
weakly affected by the periodic drive. At a critical value 
of the driving frequency (or equivalently of the driving am¬ 
plitude), the system undergoes a dynamical phase transition 
and flows towards the inhnite temperature absorbing (ergodic) 
state. The transition occurs at a finite value of the drive am¬ 
plitude and is therefore beyond the reach of perturbative ap¬ 
proaches. To provide a glimpse about the nature of the tran¬ 
sition we considered the lowest-order Magnus expansion sug¬ 
gesting that the transition could correspond to the point where 
the kinetic term in the Eloquet Hamiltonian becomes negative. 

The existence of a transition at a hnite value of the driv¬ 
ing amplitude is further supported by two numerical methods: 
a self-consistent time-dependent variational approach, and a 
semiclassical approach. Interestingly, the emergent phase dia¬ 
gram (Eig.sj^an d|3 is in qualitative agreement with a simple- 
minded quadratic expansion (Eig.[^. The dynamical phase di¬ 
agram displays two distinct stability islands, respectively for 
large driving frequencies 7 and for large go = 1/T g(t). 

The former island becomes unstable when 7 is comparable to 
the short wavelength cutoff A. This suggests that the transi¬ 
tion is of mean-held nature, thus bearing the same character as 
the parametric resonance of a single harmonic oscillator. The 
latter island is related to the existence of hnite excitation gap, 
which protects the system from low-frequency drives. In both 
cases, we observe that quantum huctuations promote ergodic- 
ity and decrease the value of the critical driving amplitude. 

The observed transitions are analogous to the stability 
threshold predicted in Ref.s II9I11I12I14I for many-body lo¬ 
calized (MBL) states, but does not require localization in real 
space. Our hndings are in contrast to the conclusions of 
D’Alessio et flZ.EH, who argued that generic many-body sys¬ 
tem should be dynamically unstable under a periodic drive. 
Ponte et al^ showed that ergodic systems are always unsta¬ 
ble to the periodic modulation of a local perturbation: we hnd 
here that global coherent perturbations can display a qualita- 
tively different behavior (See also Russomanno et al^ for a 
specihc example). 
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Appendix A: Third-order Magnus expansion 


If the driving frequency 7 is the largest scale in the problem, 
Magnus expansion applies and going up to third order one has: 


= ^£dtH{h) = 

= J dx[KP^ + - gocos{(j))] 

dh = 0 



x({F(U) , {F(/2),F(/3)}} + {F(f3),{[F(/2),F(U)}}) 


where the time-integral domain is ordered 0 </„<... < 
t 2 < ti < T and T is the period of the driving T = —. 

( 2 ) 

In the present case the second-order term vanishes = 
0 : for time-reversal invariant perturbations, all even-order 
terms are exactly zero. The third order Magnus expansion 
leads to the effective Hamiltonian (|^, where we used the fol¬ 
lowing identities: 

{P^, cos((^)}} = — 2 /{P^, Psin((/))} = 4 P^ cos{ 4 ') ; 

(Al) 


= 2{dx4>) [dx sm{(t)) P sin(^)5a:] ; 

(A2) 



[gih) - g{t 2 ) +g{ti) - 5 (^ 2 )] = 51 ; 

(A3) 



g{h)[g{h) - g{t 2 )\ + g{h)[g(ti) - g{t 2 )\ = gogi - ^gf ■ 

(A4) 


As mentioned in the text, the limit A —^ 0 of Eq. 0 re¬ 
covers the case of an isolated classical Kapitza pendulum. In 
Fig. I^we show that this limit is correctly reproduced by the 
self-consistent variational approach. With decreasing gi, the 
time required for the system to become unstable grows ap¬ 
proximately as l/gi- This explains the apparent inconsistency 
between the quadratic expansion and the self-consistent vari¬ 
ational approach (Fig.|^a)) for A /7 ^ 1 and Kgi/^'^. 



0 0.1 0.2 0.3 0.4 0.5 


Kgoh"^ 

FIG. 9: Stability diagram obtained from the self-consistent varia¬ 
tional approach for A /7 = 0.01, K = O.OStt. The color code 
shows logfyTd), where ra is defined as Z{Td)IZ{0) = 10 “^ (dark 
green areas correspond to Td > Tf — 10'^(27r/7)). Orange dashed 
lines are from the quadratic expansion and solid black lines are from 
Broer et aP^. The inset plot shows Td as a function of Kgi/~f^ 
along the resonance Kgoj^'^ = 1/4. The dashed line refers to 
Td/ ~ {Kgi/j^yK 


Appendix C: Finite time scaling 


In this Appendix we extend the analysis of Sec. |VIII and 
determine the scaling of physical observables as a function of 
the number of driving periods. Figs. [T^ a) and[TTJa) show the 
energy absorption rate and the average cosine as a function 
of the driving amplitude, at different times. These plots show 
that all the curves tend to a well defined long-time limit. To 
better highlight this asymptotic limit, we shift each curve to 
take into account the dependence of the critical amplitude on 
Tun (Fig.| 8 ). Specifically we consider finite-time corrections 
of the type: g^ —)■ gd^ + A/Tun)- Through this transforma¬ 
tion we obtain an excellent data collapse on a single universal 
asymptotic curve, as shown Fig.s [^b) and [^b). In both 
cases, the asymptotic curves become steeper and steeper as 
the number of driving periods increases. This may indicate 
that physical observables are ultimately not continuous in the 
Tun —f 00 limit, with important consequences for the univer¬ 
sal properties of the transition. 
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(a) (a) 




FIG. 10: Energy absorption rate as a function of the normalized peri- FIG. 11: Same as Fig.[T^for the amplitude of the cosine oscillations. 

odic amplitude, for a different number of periods (Tfin). (a) raw data 

for a system of size L = 400, K = O.drr and A/y = 0.1; (b) same 

data on a normalized x-axis, showing a good data collapse when the 

position of the peak is rescaled as, g = + A/T with A = 20. 
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